Percolation in networks composed of connectivity and dependency links 
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Networks composed from both connectivity and dependency links were found to be more vul- 
nerable compared to classical networks with only connectivity links. Their percolation transition 
is usually of a first order compared to the second order transition found in classical networks. We 
analytically analyze the effect of different distributions of dependencies links on the robustness of 
networks. For a random Erdos-Renyi (ER) network with average degree k that is divided into 
dependency clusters of size s, the fraction of nodes that belong to the giant component, Poo, is 
given by Poo = p s_1 [1 — exp (— fcpPoo)] s where 1 — p is the initial fraction of removed nodes. Our 
general result coincides with the known Erdos-Renyi equation for random networks for s = 1 and 
with the result of Parshani et al (PNAS, in press, 2011) for s = 2. For networks with Poissonian 
distribution of dependency links we find that Poo is given by Poo = /fc, p (Poo)e^ s>_1)(p/fc p (Poo) " i:i 
where fk, P (Poo) = 1 — exp( — fcpPoo) and (s) is the mean value of the size of dependency clusters. 
For networks with Gaussian distribution of dependency links we show how the average and width 
of the distribution affect the robustness of the networks. 



I. INTRODUCTION 

Many systems can be efficiently modeled using a net- 
work structure where the system entities are the net- 
work nodes and the relations between the entities are 
the network links [H-ll4j|. However, many systems are 
also characterized by small sub groups in which the enti- 
ties belonging to a group strongly depend on each other. 
We coin the relation between each two nodes in such a 
group as dependency links [15|]. For example consider a 
financial network: Each company has trading and sales 
connections with other companies (connectivity links). 
These connections enable the companies to interact with 
others and function together as a global financial market. 
In addition companies that belong to the same owner 
strongly depend on one another (dependency links). If 
one company fails the owner might not be able to finance 
the other companies that will fail too. Another example 
is an online social network (Facebook or Twitter): Each 
individual communicates with his friends (connectivity 
links) , thus forming a social network through which infor- 
mation and rumors can spread. However, many individ- 
uals will only participate in a social network if other in- 
dividuals with common interests also participate in that 
social network, thereby forming dependency groups. 

Previous studies focused on network models containing 
only a single type of links, either connectivity links jiol . 
Il6l - f20| or dependency links [2lT - [25j . A network model 
containing both connectivity and dependency links was 
first introduced for two interdependent networks (2(1 [27| . 

A recent paper [l5[ introduced for the first time a sin- 
gle network model containing both connectivity and de- 
pendency links. In this network model the initial failure 
of nodes may trigger an iterative process of cascading 
failures that has a devastating effect on the network sta- 
bility. The cascading failures are a result of the synergy 
between two different effects: a) A percolation process 
governed by connectivity links, b) The failure of an en- 
tire dependency groups due to a failure of one member 




FIG. 1: (Colors online) Connectivity network with depen- 
dency clusters. The edges represent connectivity relations 
while the (blue, red and pink) groups surrounded by curves 
represent dependency relations between all the nodes of the 
same group (color) . The dependency relations can be between 
very "far" nodes, in the connectivity network. In the gen- 
eral case, the sizes of each dependency clusters follow a given 
distribution. 



within the group. For a high density of dependency links 
the network disintegrates in a form of a first order phase 
transition while for a low density of dependency links the 
network disintegrates in a second order transition. 

However the combined model presented in 15] was 
based on an unrealistic assumption that all dependency 
groups are of size 2, i.e. only pair of nodes depend on 
each other. In reality, as the examples above suggest, 
groups of several elements may depend on each other. 

In this paper we analyze both analytically and numer- 
ically the general case of a network with different sizes 



2 



of dependency clusters, as illustrated in Fig. [T] We 
study networks with three different types of dependency 
groups: a) Fixed size, s, dependency groups, b) Nor- 
mally distributed sizes of dependency groups, c) Poisson 
distributed sizes of dependency groups. We find that 
for random networks with an average degree k which are 
divided into dependency groups (clusters) of size s, the 
fraction of nodes belonging to the giant component, Poo, 
is given by Poo = p s ~ l [1 — exp (— kpPooff , where 1 — p is 
the initial fraction of removed nodes. The critical thresh- 
old, p c , below which the network collapses (Poo = 0) is 
given by Eq. (fT5|) . Our result for s = 1 (a node de- 
pends only on itself) coincides with the known Erdds- 
Renyi equation, = 1 — ex p (— fcp Ppo), for a network 
without dependency relations |18l - l20l | . We also show that 
for s > 2 a process of cascading failures occurs and the 
percolation transition is of first order. 

For normally distributed dependency groups with an 
average size (s) and width a, we find that the system be- 
comes more stable (smaller p c ) for a broader size distribu- 
tion (FigJS])- When a — > the results are the same as the 
case of fixed size dependency groups with (s) = s (Figj3]). 
We also analyze both analytically and numerically the 
case of a Poisson distribution of dependency cluster sizes 
and obtain analytical equations for both Pco((s)) and p c 
(Eqs. (Ell),® and ©). 



II. GENERAL FORMALISM 

When nodes fail in a network containing both connec- 
tivity links and dependency clusters, two different pro- 
cesses occur, (i) Connectivity links connected to these 
nodes fail, causing other nodes to disconnect from the 
network (percolation step), (ii) A failing node cause 
the failure of all the other nodes of its dependency clus- 
ter, even though they are still connected via connectivity 
links (dependency step). Thus, a node that fails in the 
percolation step leads to the failure of its entire depen- 
dency cluster, which in turn leads to a new percolation 
step, which further leads to a dependency step and so on. 
Once the cascade process is triggered it will only stop if 
nodes that fail in one step do not cause additional failure 
in the next step. 

We start by presenting the formalism describing the 
iterative process of cascading failures. On each step 
we apply the two processes - a percolation process fol- 
lowed by the removal of relevant dependency groups. 
Before each percolation stage the accumulated cascades 
are described as equivalent to a single random removal, 
1 — tp%. Similarly, before each dependency stage the ac- 
cumulated cascades are equivalent to a single random re- 
moval, 1 — ipn- When applying the percolation process at 
stage n on a network of size ip? the remaining giant com- 
ponent consists of a fraction g p (ipn) which is a fraction 
'Pn — V'n5p(' ( /'n) from the original network. Similarly, ap- 
plying the dependency process at stage n on a network of 
size ip® results in a remaining functional nodes consisting 



of a fraction gD(ipn) which is a fraction (p® = ipn9D(^n) 
from the original network. 

The iterative process is initiated by the removal of 
a fraction 1 — p of the network nodes. The remaining 
part of the network is ipf = p. This initial removal will 
cause additional nodes to disconnect from the giant clus- 
ter due to the percolation process. The fraction of nodes 
that remain functional after the percolation process is 
4>i = 7 / , i5p(V'i)- Before the dependency step we describe 
the accumulated cascades of the previous steps. The frac- 
tions of nodes that fail due to the initial removal and due 
to first percolation step are 1 — ip\ and ipi — 4>\, respec- 
tively, and the accumulated cascades are equivalent to a 
single random removal of 1 — <jf[ = (1 — ipi) + {ipf — 4%)- 
Thus, we denote the remaining functional part before the 
dependency step as ip® = <p\ . Each node from the non 
functional part will cause all the other nodes of its de- 
pendency cluster to also fail (dependency process). The 
remaining functional part of the network after the depen- 
dency step is 01 = ipfgD(tfji). 

Let us now calculate the accumulated failure up to this 
step. The sum of the previous steps, the initial removal 
of (1-p), the removal due to the percolation step, (tpf — 
4>\) and the removal due to the dependency step (ip® — 
0f ), is equivalent to a single random removal of (1 — 
pg D (jPi)) from the original network ([IB]). After such 
removal the remaining part of the network before the 
second percolation step is tp^ = P9 vl ) \)i and the size of 
giant cluster is then (p? 2 = V's'SpWD- 

Following this approach we can construct the sequences 
4>n and ipn °f the remaining fraction of nodes and the 
sequences (p^ and (p® of functional nodes, at each stage 
of the cascade of failures. The general form is given by: 

< = P,<Pl = ^i^W ); 



To determine the state of the system at the end of the 
cascade process we look at ip"^ and ipm at the limit of 
to — > oo. This limit must satisfy the equation ip^ n =ip^ n+1 
(or V , m = '0m+i) since eventually the clusters stop frag- 
menting and the fractions of randomly removed nodes at 
step to and to + 1 are equal. Thus, at steady state the 
system satisfies the set of two equations 

= <to(V&)p, 

= (2) 

Denoting x = V'S and y = ip^ we arrive to a system 
of two equations with two unknowns: x — pg p (y) y — 
pgjj(x) which can be reduced to: 

x = pg P (pgD(x)). (3) 
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Solving equation (J3]) we obtain the size of the network at 
the end of a cascade initiated by random removal of 1 — p 
of the nodes. 

Next, we calculate explicitly <7d(P) and g P {T). In the 
general case, each node belongs to a dependency group 
of size s with a probability q(s) so that the number of 
groups of size s is equal to q(s)N / s. Since after random 
removal of 1 — T of the nodes each group of size s remains 
functional with a probability T s , the total number of 
nodes that remain functional is given by q(s)NT s . 

Thus, we define the function go{T) as the fraction of 
nodes that remain functional out of the TN nodes that 
were not removed, 



g D (T) = J2l( s ) TS 



(4) 



Analogous to gn(T), g P (T) is defined as the fraction of 
nodes belonging to the giant cluster of the connectivity 
network after random removal of 1 — T of the nodes. The 
percolation process can be solved analytically by using 
the apparatus of generating function. As in Refs. [28l — 
|30| we will introduce the generating function of the de- 
gree distributions Go(£) = Sfe-P(&)£ fe - Analogously, we 
will introduce the generating function of the underlining 
branching processes, Gi(£) = Gq(£)/Gq(1). Random re- 
moval of fraction 1 — T of nodes will change the degree 
distribution of the remaining nodes, so the generating 
function of the new distribution is equal to the generat- 
ing function of the original distribution with the argu- 
ment equal to 1 — T(l — £) 28]. The fraction of nodes 
that belong to the giant component after the removal of 
1 — T nodes is [H& 

g p (T) = l-G [l-T(l-u)}, (5) 

where u = u(T) satisfies the self-consistency relation 

u = G 1 [l-T(l-u)]. (6) 

III. ER NETWORKS 

The formalism presented in Sect. II is general for a ran- 
dom network having any degree distribution. In the case 
of an ER network, whose degrees are Poisson-distributed 
[l8l - l20j . the problem can be solved explicitly. Suppose 
that the average degree of the network is k. Then, 
Gi(0 = G (0 - exp[fc(£ - 1)]. Thus, g p (x) = 1 - u 
and therefore Eq. ([3]) becomes 



x =p[l -u], 
where u is defined according to ©, by 

u = e-xjp(—kpgn(x)(l — u)). 



(7) 



(8) 



Using the definition of gn(x), Eq. @, together with 
Eq. © we get the general solution for the steady state 
of the network at the end of the cascade failure process 
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FIG. 2: (a) The size of the giant cluster, (j>oo = Poo • p, vs. p, 
the fraction of nodes that remain after random removal, for 
ER networks (k = 8) for different fixed sizes of dependency 
clusters, s. The symbols represent simulation results of sys- 
tems of 50, 000 nodes and the solid lines show the theoretical 
predictions. For the case of s = 1 there are no dependency 
clusters and the regular percolation process leads to second or- 
der phase transition. For s > 2, a first order phase transition 
characterize the percolation process represented by disconti- 
nuity of Poo at p c . Both the regular and the new first order 
percolation obey Eq. (| 15f) . (b) The number of iterative fail- 
ures (NOI) sharply increases when approaching the critical 
threshold p c for the first order transitions, and thus repre- 
sent a useful method for identifying accurately the value of p c 
[l5l ]. In the figure, each curve is maximal as its related curve 
in panel (a) approaches from both sides the critical threshold. 



u = e 



-fc£r=i<?(s)p 3 (i-«) s 



(9) 



In order to present u, obtained from Eq. @, in terms 
of Pqo recall that at the steady state the size of the giant 
cluster (j)oo = (f)? = (ft®, and according to |T]) 

00 00 
$00 = xg D {x) =x^2 it^x^ 1 = X! l(s)x s - (10) 

s=l s=l 

Since = (fioo/p, we get the relation, 

00 

p o-p=Y,i( s )p s ( i - u y ( n ) 

s=l 

and using @, a simple equation for Poo is obtained 

„ In m . , 

p - - ~w (12) 

where u is the solution of Eq. ©. 
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Up to this point, we obtained the size of the network 
at each step of the cascade process, Eq. (JT|), and in par- 
ticular, its size, Poo, Eqs. ((9|) and (fl~2j) . at the end of the 
cascade, for the general case of a given distribution q(s) 
of sizes of dependency clusters. 



IV. FIXED SIZE OF DEPENDENCY CLUSTERS 

Using the general solution described above, we analyze 
the case of a fixed size, s, of dependency clusters. In 
particular, we find the size of the giant component and 
the critical fraction of the network, 1 — p c , that, if re- 
moved, leads to complete fragmentation of the network. 
In this case, go, given in (|4]), becomes 9d{T) — r s_1 
and Eqs. (j9|) and (fTTj) become respectively: 



Pn. 



P s - x (l-u) s 



which can be combined into a single equation: 

Poo = f- 1 (i - e- k p p ™y . 



(13) 
(14) 



(15) 



Eq. (Tl5|) coincides for s = 1 (a node depends only 
on itself) with the known Erdos-Renyi equation [18t - l20| . 
Poo = 1 — cxp(— kpPoo), for a network without depen- 
dency relations. Moreover, for s — 2, Eq. (IT5|) yields the 
result obtained in [l5[ for the case of dependency pairs. 

Fig. [2] shows the size of the giant cluster, P^, versus 
the fraction of nodes, p, remaining after an initial random 
removal of 1 — p, for the case of ER network with fixed 
size of dependency clusters s. The case of s = 1, each 
node depend only on itself, is the regular second order 
percolation transition. For any s > 2, a first order phase 
transition characterizes the percolation process. Both 
the regular and the new first order percolation obey Eq. 

Finding the transition point via simulations is always 
a difficult task that requires high precision. In the case 
of s > 2 where first order transition occurs we are able 
to calculate the transition point with good precision by 
identifying the special behavior characterizing the num- 
ber of iterations (NOI) in the cascading process 15;] . This 
number sharply drops as the distance from the transition 
point is increased. Thus, plotting the NOI as a function 
of p, provides a useful and precise method for identifying 
the transition point p c in the first order region. Fig. [2Jd 
presents NOI of the simulation results of Fig. [2ji. The 
transition point, p c , can easily be identified by the sharp 
peak characterizing the percolation threshold. The re- 
sults shown in Fig. 2b are in excellent agreement with 
theory. 

Next, we find analytically the percolation threshold, 
p c , for the case of a fixed size of dependency clusters. 
Eq. (JT3J) , which is the condition for a steady state, have 
a trivial solution at u = 1, corresponds, by (fl"2"|) . to a 
complete fragmentation of the network. For large p there 
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FIG. 3: Theory (Eq. (|18l) . dashed lines) and simulation re- 
sults (symbols) are compared for the values of p c for ER net- 
works with different average degree k and different fixed sizes 
of dependency clusters, s. The full symbols, for s = 1 rep- 
resent the known Erdos-Renyi second-order phase transition 
threshold for a network without dependency relations, while 
open symbols represent first-order phase transition thresh- 
olds. 



is another solution of < u < 1, corresponding to a 
finite fraction of the network. Therefore, the critical case 
corresponds to satisfying both the tangential condition 
forEq. m 



1 



(16) 



[kp s s(l 



as well as Eq. (fT3j) . Thus, combining Eqs. (|T6|) and (|T3 
we get a closed- form expression for the critical value, u c 



u c = exp 



(17) 



Once u c is known, we obtain p c by substituting it into 

Eq. eg 



P, 



= [ksu c (i - u c y 



-l/s 



(18) 



For s = 1 we obtain the known result p c = 1/fc of 
Erdos-Renyi ^M^. Substituting s = 2 in Eqs. (H7J) 
and (fl~8| one obtains p c 2 = l/[2fcu c (l — u c )], which coin- 
cides with the exact result found in [l5|. In Fig. [3] we 
plot the values of p c as a function of s for several k val- 
ues. Note the dramatic effect of the dependencies on the 
vulnerability of the system. Even for high values of k, p c 
approach rapidly to 1 even for relatively small s values. 

Next, we show that the difference between the con- 
tinuous second order percolation transition and the first 
order transition is characterized not only by the abrupt 
jump in the size of the giant cluster at the critical point, 
p c , but also by a difference in the scaling behavior of the 
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FIG. 4: The minimum averaged degree, k m in, as a function of 
the (fixed) size of dependency clusters, s. The solid line shows 
the theoretical results, obtained from Eqs. (|17[) and (|22l) . The 
region of s and k values below the line represents an unstable 
network that will collapse after any single node failure. For k 
and s values above the line, the network is stable and there 
exists p c < 1. 



giant component, </>oo, near p c . The scaling near p c is 
denned by the exponent /3, 



4>oo{p) ~ <Poo(Pc) ~ (P-Pcf- 



(19) 



Eq. (|15p can be written in terms of </>oo(= -Poo - p), the 
size of the giant component, 



cboo =p s {l-er k ^y 



(20) 



For the case of s = 1 (ER) 0oo(p) — </>oo(Pc) changes (as 
well known [13]) linearly with p — p c and (3 — 1. For 
s > 2, we calculate the behavior close to (and above) the 
critical point, 

P = p c + e 

0oo = + S 



when S, e —> and 
can be written as 



+ 5) 



l/s 



1 



where A 



-fe(0So+<5) 



sl/a 



^oo(Pc)- For this case, Eq. (|2"0")) 



= A [1 + Cii + C 2 <5 2 + . . .] 

(21) 



? iven by C x = ( -4- 



p c and the linear coefhcient is 

e fc^_i )- However, using Eqs. 
([T2"T) and (fT7|) we obtain that C\ = 0, so near the critical 
point e ~ S 2 and <5 ~ e 1 / 2 . Thus, the scaling behavior of 
the giant component near the first order transition, Eq. 
(p~9|) . is characterized by the critical exponent /3 = 1/2. 

For a fixed s, when k is smaller than a critical number 
kmin(s), p c > 1, meaning that for fc < k m i n (s), the net- 
work will collapse for any finite number of nodes failure. 
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FIG. 5: (a) The size of the giant cluster, 0oo = -Poo • p, vs. 
p (solid lines) and (j>°°(Pc) (symbols), for ER networks with 
average degree k = 15 and dependency clusters normally dis- 
tributed around averaged size (s) = 4. The different curves 
represent different standard deviations, a (a — (full triangle 
right), a — 0.8 (circle), a = 1.6 (square), a = 2.4 (triangle 
up) and a = 3.2 (diamond)), (b) Theory (solid lines) and sim- 
ulations (symbols) values of p c for ER networks with average 
degree k = 15 and dependency clusters normally distributed 
around averaged size (s), versus the width of the distribution, 
a. The three curves represent different values of js/,: (s)=4 
(circles), {s)—6 (squares) and (s)=10 (triangles). For a — > 
(full symbols) the Gaussian distribution becomes a 8 function 
with fixed size of dependency clusters and thus p c is identical 
to those obtained by Eqs. (|17|) and (|18p for groups of single 



From Eq. (|18p we get the minimum of k as a function of 

s, 



k mm {s) = [su c (l - u c ) s : ] 



(22) 



Fig. U shows the minimum averaged degree, k m i n , as 
a function of the size of dependency clusters, s. 
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V. GAUSSIAN DISTRIBUTION OF 
DEPENDENCY GROUPS 

Using the general solution, given in Eqs. (j9|) and (|TT|) . 
one can calculate after initial removal of 1 — p of 
the nodes and get p c for every distribution of sizes of 
dependency clusters, q(s). Here, we calculate p c in the 
case of a normal Gaussian distribution for the size of the 
dependency clusters with average size (s), and variance 
a 2 . In this case, the probability of a random node to 
belong to a dependency cluster of size s is given by 



Ae 



-(s-(s)) 2 /2a 2 





1 < s < 2(s) - 
elsewhere, 



(23) 




where A is a normalization constant. Note that q(s) ^ 
only for 1 < s < 2{s) — 1 in order to have a symmetrical 
distribution around (s). 

This case generalizes our results of dependency clus- 
ters (for a — > 0) having a single size s, and show the 
deviations from these results as the distribution becomes 
broader. In this case there are nodes that belong to 
dependency clusters which are larger then (s) and thus 
have higher probability to become non-functional while 
the same number of nodes belong to dependency clusters 
which are smaller then (s) and have smaller probability 
to become non- functional. We find that the first order 
transition threshold decreases as the distribution of de- 
pendency groups becomes broader and thus, the network 
is more stable, as shown in Fig. HJa). However, this effect 
becomes weaker for larger (s), as shown in Fig. [5jb) . 



VI. POISSON DISTRIBUTION OF 
DEPENDENCY GROUPS 



FIG. 6: The size of the giant cluster, 4>oo = Poo • p, vs. p, for 
ER networks (k = 8) and Poisson distribution of dependency 
clusters with different average sizes, (s). The symbols rep- 
resent simulation results of systems of 50, 000 nodes and the 
solid lines show the theoretical predictions. For (s) = 1 and 
2 the network undergoes a second order transition while for 
(s) = 2.5, 3.5 and 5 the network undergoes a first order tran- 
sition (see figure (J7J where the exact transition point from 
first to second order transition is shown). 



A 

i — 

n 



lim 



A(T-1) 



e A(T-l) = e (( s )_l)(T-l)_ 



Following (|7|) and © , we get an equation for u for the 
case of Poisson distribution: 



Next, we study the case of a Poisson distribution of 
dependency cluster sizes. In this case, the probability 
that a random node is depended on s' other nodes (and 
thus the size of the cluster, s, equal to s' + 1) is 



e 



n->oo \s' I \n 



A 



(24) 



when A = (s) — 1 is the average number of other nodes 
that depend on a random node. The dependency process 
can be calculated, using Eq. (0} and considering that 
q(s) =p{s - I), 



9d 



(T) =^Tp(s- 1)T S_1 = £ p{s')T s \ (25) 



s=l 

with s' = s — 1. Usine 



s'=0 

we obtain, 



9d(T) = £ lim ( n ) ( * 

* — ' n-s-oo \s' I \ n 



'=0 



n 



bxu = -kp(l - uJefW-ilW 1 -")-!). (26) 
Finally, using (IT21 is obtained in a closed form, 



«^)-l)(p/ fc ,p(P=o)-l) 



(27) 



where fk, P { p ^>) = 1 _ ex P (~kpPoo)- 

Fig. IH1 shows the size of the giant cluster at steady 
state versus p for different values of (s). For small (s), 
many nodes do not depend on other nodes so the effect 
of the dependency clusters is rather weak and thus the 
percolation transition is of second order, meaning that 
Poo continuously decreases from a finite value to zero at 
a specific transition point, p 1 * . However, for large (s) 
the network undergoes a first order transition, meaning 
that the size of the giant cluster abruptly jumps discon- 
tinuously from finite size for p > p\ to zero for p < p\. 
Such a network is qualitatively more vulnerable than a 
network that undergoes a second order transition due to 
the cascading failures leading to high vulnerability of the 
network around p\. 
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which means that the network is completely fragmented. 
The second order transition point, p l * , corresponds to 
the solution of Eq. ([25)1 where u —> 1. This condition 



gives p, 



ii 



P. 



ii 



(28) 



Note, that for the case of (s) = 1, meaning that all the 
nodes are not depended, p 1 * — 1/fc as for regular Erdos- 
Renyi network. 

The first order transition point, p^, corresponds to the 
tangential intersection of the left and right terms of Eq. 
(121)1) , meaning that the derivatives of both with respect 
to u are equal. This yields, 
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FIG. 7: (a) Theory (lines) and simulations (symbols) are com- 
pared for the values of pl({s}) and pi 1 ((«)) for ER networks 
with different average degree, k (circles for k — 3, triangles 
for k — 5 and squares for k = 8). For (s) > (s) c - (dashed 
line) the network undergoes a first order transition. The the- 
oretical values of the transition point, pl{(s}), that are cal- 
culated according to Eq. I|29p are compared with simulations 
(open symbols) performed using the NOI method (explained 
in text). For (s) < (s) c the network undergoes a second order 
transition, and the theoretical values of the transition point, 
p I c I ((s}), that are calculated according to Eq. (|28p are com- 
pared with simulations (full symbols) performed using the 
second largest cluster method. The dashed line separating 
between the first and second order is obtained according to 
Eq. (|32[) . (b) The size of the giant cluster, above the crit- 
ical point, p c , is described by (/>oc(p) ~ (pcc(Pc) ~ (p-p c Y , as 
shown in the inset. Critical exponent, /3, above the transition 
point, Pc, is plotted versus the average size of the dependency 
clusters, (s), for a network with k — 8 (compare to panel 
(a); squares). For the region of the second order transition, 
(s) < (s)c, we find /3 = 1 ((s) c is marked by the vertical 
dashed line) while for the region of the first order transition, 
(a) > (s) c , P = 1/2. 



Next, we find explicitly, by analyzing Eq. (I26I) . the first 
order transition point, and the second order transition 
point, pi 1 . Eq. (|2l)]) has a trivial solution for it = 1, 



u - 1 



j In '( 



(29) 



where u is the solution of Eq. (|26| . 

Fig. [7J (a) shows p\ and p 1 / versus the average size of 
the dependency clusters, (s). At critical values p = p* 
and (s) — (s) c the phase transition changes from first or- 
der to a second order. The values of p* and (s) c are ob- 
tained when the conditions for both the first and second 
order transitions are satisfied simultaneously. Applying 
both conditions we obtain 



2«s)c-l) 



1 

P*c 



(30) 
(31) 



For a given ER network with average degree fc, Eq. 
([3TJ| provides the critical average size of dependency clus- 
ters, (s) c . A network with (s) < (s) c undergoes a sec- 
ond order phase transition while for (s) > (s) c the net- 
work undergoes a first order transition. Therefore, p*, 
obtained from Eq. (1311) for the case of (s) = (s) c , charac- 
terizes the stability of a network with maximal (s) under 
the constraint of undergoing second order transition. The 
critical case, shown in Fig. [7] (a), of transition from first 
order to second order transition (dashed line) is given by 



1 



Pc 



2«*>c-l)" 



(32) 



As shown in Figs. [S] and [7J increasing the size of the 
dependency clusters increases the network vulnerability 
and p^ becomes larger. A critical average size of depen- 
dency clusters, (s) max , corresponds to p J c = 1 meaning 
that the network completely fragments as a result of re- 
moval of even a single node. Such a network can be 
regarded as unstable. The value of (s) max is given by 
Eqs. (j29|) and (l26l) under the condition of pi = 1. 

Thus, the stability of a random network with depen- 
dency clusters with average size (s) with Poissonian dis- 
tribution can be summarized 
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(s) < (s) c Second order transition 

(s) c < (s) < (s) max First order transition 

(s) > {s)max Unstable network. 
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